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A definition of an effective meson-meson interaction adapted to the framework of a lattice simulation is pre- 
sented. Results, based on a truncated momentum-space 4-point time correlation matrix, and preliminary data 
from a complementary coordinate-space simulation, are shown. 



1. INTRODUCTION 

The nucleon-nucleon interaction has been 
parametrized by boson-exchange models for se- 
veral decades. By their very ansatz, boson- 
exchange models lack the possibility of giving in- 
sight into the nature of the nucleon-nucleon inter- 
action and its underlying physics from first prin- 
ciples. 

From the viewpoint of lattice QCD a study 
of effective hadron-hadron interaction opens up 
mainly uncharted territory. Only a few attempts 
have been made to break new ground [1,2] in 
a field which could be characterized as Lattice 
Nuclear Physics. This field would concern itself 
mostly with the effective forces between composite 
colour singlets, in other words, their "chemistry". 
This goes along with a characteristic change of 
scale leaving a residual effective interaction 10 -2 
to 10 -3 times smaller than a typical hadron mass. 

The numerical difficulties involved in going 
from the GeV to the MeV scale on the lattice are 
evident. For this reason it seems appropriate to 
investigate simpler lattice models which exhibit 
similar physics. Most importantly, of course, one 
needs to insist on confinement and the presence 
of fermions. 

It is in this spirit that we report on a simu- 
lation aimed at extracting an effective meson- 
meson interaction from a QED gauge field model 
in 2 + 1 dimensions with staggered fermions which 
has confinement. In an earlier related simulation 
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[3] scattering phase shifts were obtained for this 
system based on Liischer's proposal. Regarding 
the parameters of the QED2+1 model and numer- 
ical strategies we refer the reader to Ref. [3]. 

2. MESON-MESON CORRELATOR 

An appropriate set of operators can be con- 
structed from one-meson fields with momentum 
P= ^(*i,*2) 

4> f (t) = L- 2 J2 e l?i x d {£t)xu(£t) (1) 

X 

where x an d X are staggered Grassmann fields 
with external flavours u and d. Meson-meson 
fields with total momentum P = then are 

= 4>-p(t)4>+p(t) • (2) 

The information about the dynamics of the 
meson-meson system is contained in the 4-point 
time correlation matrix 

c$(t,t ) = <4«<MM> - KMX^Co)) • (3) 

Working out the contractions between the Grass- 
mann fields the following diagrammatic represen- 
tation is obtained 

C (4) _ C (4A) + C (4B) _ C (4C) _ C (4D) ^ 

= 1 1 + X - M - M • < 5 > 

The quark-exchange diagrams C^ AC ) and C^ AD ) 
are a source of interaction between the mesons, 
but then, so are diagrams C* 1 - 4 ^- 1 and C^ AB \ De- 
noting contractions of the Grassmann fields by 

n n 

X...x=G n with n = l..A (6) 
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we have for example 

43 21 12 34 

c (4A) = ( 4 ^ V-^> ~ (G1G2G3G4} (?) 

where ~ stands for four Fourier sums that carry 
over from (1). Note that the gauge configuration 
average is taken over the product of all four prop- 
agators G. This is a source of gluonic interaction. 

In order to define an effective meson-meson in- 
teraction it is crucial that the noninteracting com- 
ponents in G 1 - 4 - 1 be isolated. This is achieved by 
means of a cumulant (or cluster) expansion [4] of 
the gauge field average. For example 

G^ ~ (GiG5)(G 3 G5) + (GiG 3 )(G5G5) (8) 

+{GiG* A ){GsG* 2 ) + (GiG^GsG^c 

where the last term defines the cumulant. The 
first term is readily interpreted to mean free, un- 
corrected, mesons on the lattice. With reference 
to (6)-(8) define 

21 12 43 34 
<? (4A) = (<f>l^_ f )( 4 <f> f ) ~ {GiG* 2 ){G 3 G* 4 } . (9) 

Diagram 

C (4B) 

can be analyzed in a similar fash- 
ion, which then leads us to define C^ AB \ The sum 
of those 

C(4) = c^ A ) + G( 4B ) (10) 

constitutes the free meson-meson correlator with 
the relative interaction switched off. 

3. EFFECTIVE INTERACTION 

The deviation of G 1 - 4 - 1 from G 1 - 4 - 1 somehow con- 
tains the residual effective meson-meson interac- 
tion. To find its definition, one may calculate 
the correlation matrix G for an elementary boson 
field <j){xt) that lives on the lattice sites and is 
subject to a perturbation, H = Ho + Hi. It is 
possible to derive an explicit expression for Hi in 
terms of the correlators of order N = and N = 1 
from the perturbative expansion of G. This result 
leads us, by way of analogy, to define an effective 
meson-meson interaction as 

%I = _iL \c^-h c ^c^] . (ii) 

dt l lt=t v ' 

This definition is consistent with C = e -^i(*-*o) 
where C is the matrix inside [. . .] in (11). 



4. RESULTS 

The simulation has been done on an L 2 x T = 
24 2 x 32 lattice with the compact U(l) Wil- 
son action at [3 = 1.5 and quenched staggered 
fermions. The correlation matrices G 1 - 4 - 1 and G 1 - 4 - 1 
were computed for a truncated set of momenta 
p = ^(^1,^2) with |fci 2| < 2, and for time slices 
t = 15 ... 19 around the symmetry point t c = 17. 
Invariance of the correlators G 1 - 4 - 1 and G 1 - 4 - 1 under 
the lattice symmetry group 0(2, Z) was utilized. 
Its irreducible representations T = A\, A2, B\, B2 
(see [3]) render the reduced matrix G 1 - 4 ' 1 "- 1 diag- 
onal, and positive, which makes computing the 
inverse square root C^ A ' V )~~ trivial. The time 
derivative in (11) was computed from cosh-fits to 
the eigenvalues of Cp 4 q F \t, to). Finally, a Fourier 
transformation, back to coordinate space, and a 
partial-wave projection were performed 
00 

{r\ni\r< ) = $>f(r, r')P t ^ -?') + ••• (12) 

1=0 

where Pt(<p) oc cos{£<p) and r, <p are the polar co- 
ordinates of r. The s-wave effective interaction 
(£ = 0) is the dominant contribution to (12). Its 
elements Ti^\r,r') are displayed in Fig. 1. Re- 
pulsion at small relative separation and attraction 
at larger distance are evident. This is consistent 
with the results of [3]. It should be noted that 
the small momenta considered in the correlation 
matrices cannot resolve short distances well. In 
particular, the repulsive core inevitably comes out 
soft and broad. 

5. COORDINATE-SPACE HARD CORE 

A complementary simulation which starts di- 
rectly from coordinate-space matrix elements of 
the correlators G 1 - 4 - 1 and G 1 - 4 - 1 reveals an inter- 
esting perspective on the repulsive core. Those 
correlator matrix elements are constructed from 
meson-meson fields with relative distance r 

Mt) = Y,e- if U-?(t)M t )- (!3) 

p 

From the resulting expression for G 1 - 4 - 1 in terms of 
the fermion propagator it is easy to see analyti- 
cally that Cii](t,t ) = if f = or f' = 0. To 
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Figure 1. Nonlocal s-wave effective interaction 
between two mesons. The vertical scale is in units 
of 10" 3 a- 3 . 



the extent that the energy W(r) of the meson- 
meson system at fixed relative distance f can be 
extracted from the large-t behavior of the diago- 
nal elements 

Cg(t,t )~c(r)e- w ™-t°) (14) 

we then may conclude that W(f = 0) = +oo, 
provided the strength factor c(r ) is non-zero at 
r = 0. Thus the effective interaction possesses a 
hard repulsive core. It is due to the anticommut- 
ing nature of the fermion field (Pauli repulsion). 
In Fig. 2 we show raw simulation data for W(r). 
A careful interpretation of those results will be 
forthcoming [5]. It should be emphasized that 
Fig. 1 and Fig. 2 are based on incompatible ap- 
proximation schemes and cannot be compared in 
quantitative terms. 

6. CONCLUSION 

A practical method to extract an effective 
hadron-hadron interaction has been developed 




Figure 2. Fixed-distance meson-meson energies 
including all four diagrams of (5). The dashed 
line indicates the mass of two free mesons. 



and applied within the framework of a QED2+1 
lattice model. (It replaces an earlier proposal [6] 
which has proved numerically unsuitable.) Its 
momentum-space version yields a fully dynami- 
cal effective interaction with short-range repul- 
sion and intermediate-range attraction, the hall- 
mark features of the nucleon-nucleon interaction. 
A complementary coordinate-space approach in- 
dicates a hard repulsive core for this model. 
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